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Abstract 

The new data on pion production in p+C interactions from the NA49 experiment at the 
CERN SPS are used for a detailed study of hadronization in the collision of protons with 
light nuclei. The comparison to the extensive set of data from p+p collisions obtained with 
the same detector allows for the separation and extraction of the projectile and taiget con- 
tributions to the pion yield both in longitudinal and in transverse momentum. 
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1 Introduction 



Hadron-nucleus interactions constitute a unique laboratory for the study of hadronization 
in the non-perturbative sector of QCD by offering, in particular, access to multiple hadronic 
collision processes and by opening an important link to nucleus -nucleus reactions. In fact the 
latter process cannot be understood properly without reference to the more elementary hadron- 
hadron and hadron-nucleus interactions. Past attempts in this field have been largely frustrated 
by the absence of detailed and consistent data sets both from hadron-hadron and from hadron- 
nucleus collisions. The NA49 experiment aims at providing such data sets obtained with the 
same detector over a wide range of projectile and target combinations [IJ by measuring yields 
of identified secondaries with a dense coverage of the available phase space. 

A recent publication concerning pion production in p-i-p collisions fT\ has now been fol- 
lowed by data of similar quality from minimum-bias p-i-C interactions [ 3 1 . This allows for a new 
attempt at comparing both reactions and extracting quantitative information about hadroniza- 
tion involving light nuclei. The hadronization process in hadron-nucleus interactions consists 
of three distinct components: 

- The projectile fragmentation which carries the imprint of multiple collisions in nuclear 
matter and is the part of principle interest for this study. 

- The fragmentation of the participating target nucleons. 

- The contribution from intranuclear cascading which is generated by the propagation and 
interaction of the participating nucleons and produced hadrons inside the nucleus. 

The separation of the three components is non-trivial if it is to be conducted in a model- 
independent way. It has been demonstrated by NA49 for the first time in p-i-Pb collisions con- 
cerning net proton production [1 1. In this case the full range of projectile and target combina- 
tions available to the experiment has been used, with baryon number conservation as the only 
external constraint. 

For pions the situation is more complex as there is no conservation law (except charge 
conservation) to be exploited. In addition the target cascading part extends closer to the cen- 
tral region than for protons, up to a; values of about -0.05. Both the target and the projectile 
components feed-over into the respective opposite hemispheres with a priori unknown ranges 
and shapes of the longitudinal momentum distributions. These parameters should be determined 
experimentally without taking recourse to unfounded model assumptions. 

The present study will attempt the quantification and separation of the relevant compo- 
nents, with a special view to the extraction of the projectile contribution in its longitudinal and 
transverse momentum behaviour. For this aim the reference to the elementary hadron-proton 
interactions is absolutely essential and most of the argumentation will be established in direct 
comparison to these processes, invoking isospin symmetry as an important constraint for the 
constitution of the target contribution. 

In its first part the paper deals with px integrated xp distributions since the superposition 
process is to first order a matter of longitudinal momentum space. After a presentation of the 
longitudinal pion yields in relation to p-i-p collisions in Sect.|2l constraints on the overlap of tar- 
get and projectile components from elementary interactions are developed in Sect.|3J resulting in 
a two-component picture of these collisions. This two-component picture is applied to the p-i-C 
reaction in Sect. H] which provides detailed information about the target and projectile fragmen- 
tation. Nuclear aspects are discussed in Sect. |5l followed by the quantification of the cascading 
contribution in Sect. |6l This allows for the extraction of the multiple collision part of the inter- 
action for a number of typical observables, and a comparison to the centrality-constrained data 
available from the experiment in Sect.|7] 
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In its second part the paper will exploit the more detailed behaviour of double inclu- 
sive cross sections which give access to the study of transverse momentum behaviour. Sect. [H 
presents the general px dependencies in relation to p+p interactions and Sect. |9l addresses in 
particular the behaviour towards higher up to about 1.8 GeV/c. A discussion of the mass and 
Pt dependence of the feed-over mechanism in relation to resonance decay will close the paper 
in Sect. 10. 



Parti 

Pt Integrated xf Distributions 



2 Comparison to p+p interactions 

2.1 Charge-average xp distributions 

The most general reduction of the multiparticle phase space of pion production that may 
be obtained is the px integrated, charge- averaged ((vr) = (tt"*" -\-'k~)/2) longitudinal momentum 
distribution, in the following characterized by the scaling variable 



Xp 



PL 



(1) 



where pl and a/s are given in the nucleon-nucleon cms. This distribution is compared to the 
p-i-p interaction |l2| by forming the ratio 



R{xf) 



{dn/dxp) 



pC 
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and shown in Fig.[T] Here the integrated pion densities dn/dxp are obtained from the interpo 
lated double differential cross sections as described in detail in [|2J and tSJ. 
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Figure 1: Ratio -R(xi;') as a function ofxp for the charge- average pion (vr) 
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Since the acceptance for particle identification of the NA49 detector imposes a limit Q at 
xf = -0.03 for the complete integration of the pt distributions, a slight extrapolation to a;^ = -0. 1 
is also indicated in Fig.[T] This extrapolation uses an extension of the data interpolation devel- 
oped for thep-i-C data [3| into the non-measured region ofpx- It will be justified in the following 
sections of this paper. The upper limit of the comparison at xp = 0.5 is due to the limited event 
statistics in the p-i-C data sample. 

Several features of interest emerge from this first comparison. The ratio R decreases 
smoothly from about 1.6 at the limit of xp = -0.1 via 1.38 at Xi? = to a value of 0.85 at 
xp = 0.5. The factor 1.6 is close to the expected mean number of projectile collisions (u) in- 
side the Carbon nucleus (see Sect. |5l) and the value 1.38 at = is close to the expectation 
((z/) -I- l)/2 from the most straight-forward projectile-target superposition scheme [4J. The de- 
crease of particle density in the far forward direction recalls a known phenomenon of projectile 
fragmentation. It corresponds to the xp dependence of the exponent a in the parametrization 

^A- (3) 



dxp 

where A is the nuclear weight. For a limited range in transverse momentum and for the far 
forward direction alpha has been shown to decrease as a function of xp 10. 

The increased particle density in the region -0.1 < xp < 0.2 on the other hand is due to the 
superposition of the two principle components of the hadronization process, namely the target 
nucleons participating in the collision and the projectile itself. The quantification and extraction 
of these two basic ingredients to the p-i-A interaction is the main aim of this study which can 
for the first time make full use of the important transition region between the two components. 
This transition region will be shown to be located between the limits -0.1 < xp < 0.1. 

A first indication concerning the validity of the extrapolation towards xp = -0.1 used in 
Fig. d is obtained by studying the sensitivity of the integrated yields to a lower cut-off in pp as 
shown in Fig. El 
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Figure 2: Ratio R{xp) as a function of xp for the charge- average pion (vr) obtained with lower 
integration limits at pp = 0, 0.2 and 0.4 GeV/c 
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Here the ratio R{xf) is obtained for lower integration limits at px = 0, 0.2 and 0.4 GeV/c. 
For the latter value the full region down to xj? = -0.1 is experimentally available. It appears 
that the effect of the successively increasing cut-offs is on the few percent level for negative 
xp and can therefore be considered as a second order phenomenon. It finds its explanation 
in the pr dependence of the range in Xi? of the projectile feed-over as discussed in Sect. |9l 
The quantitatively larger effect in the forward hemisphere is noteworthy. It has its origin in the 
sizeable increase in (px) for p-i-C compared to p-i-p interactions with ultimately very important 
Pt dependencies in the range pr > 1-5 GeV/c (Sect.lU). 

2.2 Charge dependence 

The yield ratio R{x p) is plotted separately for 7r+ and tt" secondaries in Fig.|3] Evidently 
there are differences both in the far forward, the projectile hemisphere, and more pronounced, 
in the backward region. The former effect results in a sizeable increase of the tt+Ztt^ ratio for 
xp > 0.3 and signals an important modification of the projectile fragmentation mechanism. The 
latter effect is directly related to the target. Due to isospin symmetry [6| the neutron content of 
the Carbon nucleus enhances the tt^ and suppresses the 7r+ yields from the target relative to the 
p-i-p interaction. The observation of the tt+Ztt" ratio in this region can therefore provide direct, 
model-independent access to the relative target contribution. 




Figure 3: Ratio R{xp) as a function of xp without pp cut off (two upper panels) and with pp 
cut off included (two lower panels) for a) and c) vr^, b) and d) vr" 
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The ratio R for different lower cut-offs in the pr integration is also shown in Fig. |3l For 
7r+ a local maximum of R at Si? ~ -0.05 is followed by a decrease towards lower x f ■ For vr^ on 
the other hand a steady increase of R through the same region of Xi? is observed. These effects 
will be studied quantitatively in Sect.|4| 



2.3 Charge ratios 

The effects discussed above result in a strong backward-forward asymmetry of the charge 

ratio 

Rc{xf) = j^^-ji — — r (4) 

for the p-i-C interaction as compared to p-i-p. This is demonstrated in Fig. H] which shows Re 
for both the p-i-C and the p-i-p reaction as a function of xf- An enlarged view of the central 
and backward xf regions is also presented in Fig.|4|for different px cut-offs in the transverse 
momentum integration, giving further support to the extrapolation used at low px- 




Figure 4: px integrated tt+Ztt ratio as a function of in p-i-p and p-i-C a) without pt cut off and 
b) with Pt cut off included as the data points in p-i-p collisions are plotted only up to = 0.02 



Again invoking isospin symmetry, an approach of Rc to unity will indicate the region of 
prevailing target fragmentation, whereas an approach to the values measured in p-i-p collisions 
will signal the prevalence of projectile fragmentation. Indeed the plots of Fig. HI allow to argue 
that both expectations are fulfilled in the regions xf < -0.1 and xf > -i-O.l, respectively. The 
situation is, however, more complex in detail. The approach to i?c = 1 is not smooth but shows 
structure with a local maximum at Xi? ~ -0.05, and the Rc values for positive xf approach but 
do not quite reach the ones measured in p-i-p collisions below xf = 0.2. For xf > 0.3 on the 
other hand, they start to exceed these values. A more detailed discussion of these second-order 
effects will be given in Sect.lHbelow. 
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3 Projectile and target contributions in elementary interactions 

The superposition and separation of the target and projectile components in p+A col- 
lisions is a main ingredient for the argumentation of this paper. The inspection of these two 
components also in the elementary hadron-proton interactions is therefore a necessary condi- 
tion for the model-independent interpretation of the data. The discussion will be built around 
three experimental facts: 

- The absence of long-range two-particle correlations at \xf\ > 0.2. 

- The presence of forward-backward multiplicity correlations at Ixj?! < 0.1. 

- The xp dependence of the tt+Ztt^ ratio in (7r)+p collisions. 



3.1 Long-range two particle correlations 

A high precision study of the forward-backward correlation function 



Q{xi,X2) 



g{xi,X2) 



-atr 



(5) 



f{Xl) ■ f{x2) 

has been performed at the CERN ISR for different combinations of identified secondary 
hadrons in p-i-p interactions. Here xi and X2 are the Feynman x values of the two particles 
which are selected in the opposite hemispheres of the reaction and / and g denote the invariant 
one and two-particle cross sections, respectively. Typical plots of Q{xi,X2) are shown in Fig.|5l 
for pion pairs and for Kvr and pvr combinations. 



O 



1 

0.5 



O 



1 

0.5 



O 



0.5 



TC+ 71"" 0.25 < < 0.3 



b) 



-T>- 



71' 71+ 0.25 <x^(^-) < 0.3 



c) 



-| T_J__L__. 



71" 71' 0.25 < < 0.3 



1.5 



0.5 



-i — 1 — , — 1 — I — , — 1 — 1 — |- 



d) 



I 

V 



K"" 71+ 0.5 < x^. < 0.55 

I 1 1 I 1 , 1 I , 1 1 I 1 L. 



0.2 0.4 0.6 0.8 1 
1.5 



1 



0.5 



e) 



^ T f T ^ 



7t+ p 0.2 < x^. < 0.25 



0.2 0.4 0.6 0.8 1 0.4 0.6 0.8 1 



Figure 5: Forward-backward correlations in p+p collisions measured by Q for: a) vr+vr"'", b) 
7r^7r+, c) 7r~7r~, d) K+7r+, and e) vr+p 
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The acceptance and particle identification constraints for the two spectrometer arms of 
the ISR experiment limit the longitudinal coverage to \xf\ > 0.2 for pions and la;^! > 0.4 for 
protons and kaons. Within these limits, no deviation of Q{xi, X2) from unity is visible within 
tight error bars with the only exception of the regions of large x p where energy-momentum 
conservation will finally impose some correlations of kinematic origin [7|. 

Some of the important results of this experiment concern the absence of charge or flavour 
exchange, the indication of gluon exchange as the exclusive exchange mechanism in soft 
hadronic collisions at the SPS/ISR energy, and the proof that no notable forward-backward 
correlations exist above the limit of = 0.2. 

3.2 Forward-backward multiplicity correlations 

A series of bubble chamber experiments |l8l|9l[T0| at SPS energies, and a streamer cham- 
ber experiment at the ISR [ 1 1 1 have extended the study of |^ to full phase space by establish- 
ing multiplicity correlations between the forward and the backward hemisphere. In experiments 
[8-10] the complete xp region is accessible by momentum measurement in the fixed-target lab- 
oratory frame. Since there is no particle identification available, the assumption of pion mass 
for the necessary Lorentz transformation creates a certain but practically negligible bias in the 
determination of central rapidity. Experiment Hill has no momentum measurement and only 
determines pseudo-rapidity r] but the 1] = reference is uniquely given for all particles by the 
angle measurement alone due to the collider geometry of the ISR. 

Forward-backward multiplicity correlations are quantified by splitting the charged multi- 
plicity rich per event into a forward (n j) and backward (rih) sample around y = or r/ = and by 
plotting the average backward multiplicity (n^,) as a function of as shown in Fig.|6l 

Except for some odd-even fluctuations at small multiplicities, a linear correlation of the 

type 

{ub) = a + h-Uf (6) 

emerges for all experiments, with values of the slope parameter h between 0.17 and 0.22 in 
the -y/s range from 14 to 24 GeV (Fig.|6^-d). These results do not contradict the absence of 
correlations shown by |7| above. In fact the observed multiplicity correlations vanish if only 
particles above one unit of rapidity or pseudo-rapidity are allowed to contribute to the sample, 
as presented in Fig. |6f-h. For the region of i/i covered and for transverse momenta in the 
neighbourhood of {pp) this corresponds to a cut mxp between 0.04 and 0.06 units as shown in 
Fig.|6^. This limit is well below the xp acceptance of \J}. 

The combined results of Sects. 13.11 and l3^ allow the statement that in the SPS energy 
range there is a limited feed-over of (mostly) pion production from the forward to the back- 
ward hemisphere and vice-versa. This feed-over contains between 17 and 21% of the particles 
produced in each hemisphere and is limited to the region of \xf\ ^ 0.05. There is no further 
long-range correlation beyond this limit. This statement has to be seen within the experimental 
uncertainties of Fig.|6f,g, and h. A tail of feed-over up to Xi? ~ 0.1 as it is following from the 
study of charge ratios. Sect. 13. 31 will correspond to a percent effect on the slope parameter h and 
is compatible with the error bars. 

Before coming to a further quantification of this situation it is worth while commenting 
on the energy dependence of the correlation determined by [l l j. In fact the correlation factor 
b increases steadily with i/i and reaches a value of 6 = 0.31 at i/i = 62 GeV. It has been 
established ifTTl that this increase is not accompanied by an extension of the feed-over range in 
xphy showing that the correlation vanishes by replacing the cut at jr^l > 1 by an s-dependent 
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Figure 6: Forward-backward multiplicity correlations, a) to d) experimental data and fits, [8- 
11]; f) to h) correlation with rapidity selection as given, [9-11]; d) xp as a function of pt at 
?/ = 1 for beam momentum between 205 and 306 GeV/c 



one at|?7|>0.51ns — 2. This cut happens to correspond in first order to an s-independent xp 
cut at 0.05. The increase in the correlation is therefore due to the non-scaling increase of central 
pion production characterized by the rising "rapidity plateau" 1 12|. 

3.3 tt+Ztt ratio in pion-induced interactions 

The range and shape of the feed-over components of hadronic interactions can also be 
constrained by using asymmetric collisions and by exploiting a measured quantity that exhibits 
a predicted behaviour for either the projectile or the target fragmentation. 
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Such a combination is given by using the tt+Ztt^ ratio for (7r)+p interactions. Here (7r)+p 
indicates the average of 7r++p and vr^+p collisions. In this case, the tt+Ztt^ ratio is exactly unity 
for the complete projectile fragmentation. This follows in a completely model-independent way 
from the fulfillment of isospin symmetry in hadronic collisions. By plotting the tt+Ztt^ ratio as a 
function of xp one therefore obtains directly the extent and shape of the proton target feed-over 
into the projectile hemisphere using the deviation of the measured charge ratio from unity. 

Experimentally this means the performance of precision measurements of identified pion 
yields in the central region both with 7r+ and vr^ beams. This has been achieved with the NA49 
detector [1 1. It has been shown that isospin symmetry, i.e. the inversion of 7r+ and vr^ yields 
when passing from a 7r+ to a tt^ projectile, is indeed fulfilled for Xi? > 0.1 in this reaction [6|. 

In the determination of the tt+Ztt^ ratio all experimental corrections drop with the excep- 
tion of the feed-down from weak decays of strange particles for the proton target contribution. 
This feed-down which is in principle asymmetric in its 7r+ and tt^ yields has been quantified. It 
contributes less than 0.5% to the tt+Ztt^ ratio in the forward region and has been corrected for. 

The resulting pionic charge ratio is shown in Fig.|7]as a function ofxp for (7r)-i-p collisions 
in comparison to the same quantity in p-i-p interactions. 
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Figure 7: tt+Ztt^ ratio in p-i-p and (7r)-i-p 

For (7r)-i-p the tt+Ztt^ ratio levels off to unity at xp above 0.08 and stays constant at this 
value up to the limit of data extraction at xp = 0.3. The same ratio for p+p collisions shows 
the known strong increase with xp which is reflecting the presence of the proton projectile. 
Surprisingly, the approach of the charge ratio to unity is not smooth but shows a characteristic 
structure with a secondary maximum at Xi? ~ 0.04 which is reminiscent of the similar behaviour 
of this ratio in the target hemisphere of p+C interactions, Fig.|4| It should be emphasized here 
that the precision of the charge ratio as far as statistical and systematic errors are concerned, is 
on the sub-percent level in this region. 
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3.4 A two-component picture of pion production in p+p and (7r)+p collisions 



The combined experimental evidence collected in the preceding sections allows for the 
construction of a two-component picture for pion production in elementary hadronic interac- 
tions. The two components correspond to the target and the projectile contributions. These con- 
tributions obey factorization, a well known phenomenon in hadronic reactions which specifies 
that the target component is independent of the projectile particle type, and vice-versa for the 
projectile fragmentation. Both components feed-over into the opposite hemisphere with a range 
that has a maximum extent of 0.1 units of in the SPS/ISR energy range. Sect. 13.31 Fixing this 
range, using the symmetry constraint at = in p-i-p collisions and limiting the integral of par- 
ticle yield in the opposite hemisphere to about 19% of the total yield per component. Sect. 13.21 
there is very limited freedom in the construction of this two-component picture. Fig. [8] shows 
the xp dependence of the two components as it is obtained obeying the above constraints, and 
their superposition to the measured total inclusive yield for p+p collisions. 



X 
■D 




Figure 8: Two-component picture of charge-averaged pion production in p-i-p collisions, show- 
ing the symmetric contribution from the target and projectile and their sum corresponding to 
the data 



The behaviour in the feed-over region -0.1 < xp < 0.1 can be further quantified by 
plotting the ratio of the proton target component relative to the total inclusive pion yield 



(dn/dxp) 



two-comp 



two-comp 



(7) 



{dn/dXF)mc\us 

as presented in Fig. 13 

In this Fig. two slightly different shapes are indicated. They correspond to the two 
pion charges. This shape difference is necessitated by the structure in the ti~^/ti~ ratio around 
xp = 0.04 observed in Fig. |71 In fact there is no a priori knowledge of the behaviour of the 
two charge states in the feed-over region: they may and indeed do differ from each other. This 
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difference is quantified in Fig.|5l) and c with respect to the charge- average distribution. It cor- 
responds to a decrease followed by a steep increase of the tt+Ztt^ ratio as one moves from the 
target hemisphere into the projectile region. This is shown in Fig. [10| which gives the corre- 
sponding prediction for the measured tt+Ztt" ratio both for p+p and (7r)-i-p collisions. 

It is interesting to note that the measurement for both reactions is precisely described, 




Figure 10: tt+Ztt ratio in p-i-p and (vr)-i-p with corresponding predictions for the proton target 
and proton projectile components 
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yielding a simultaneous explanation for the slight shoulder and the structure around xp = 0.04 
observed in p+p and (7r)+p collisions, respectively. This corresponds to a long-range charge 
correlation within the target and projectile contributions with the respective incoming hadron, 
which is for the first time extracted here. This correlation will present an interesting challenge 
for microscopic production models, in particular for those based on string fragmentation ideas, 
as the memory of the original parton charge tends to be quickly lost along the string lfT3ll . 

3.5 Comparison to the feed-over in net baryon production 

The NA49 data from p+p and (7r)+p interactions have been previously used to extract 
the two-component nature of net baryon production. The method used is different from the 
one applied for the pion production as net baryon number conservation can be invoked as a 
powerful constraint. The argumentation gives consistent results using pion induced collisions 
and the forward-backward baryon correlation in p-i-p events [14|. The extracted shape of the 
feed-over distribution is directly and independently measured by three methods in a completely 
model-independent way. The result is shown in Fig. [TT]in comparison to the pion feed-over 
obtained above. 




Figure 1 1 : a) Target component for net protons measured with different methods and b) target 
components of (tt) and net protons 

There is clearly an important difference between net protons and pions concerning the 
range of the feed-over. This indicates a strong mass dependence. In connection with the yield 
increase with transverse momentum, Sect.[8l it will be argued that there is also a p-r dependence 
of the feed-over. Both the mass and the pt dependence will be shown in Sect. 10 to follow from 
resonance decay. 

4 Projectile and target components in the p+C interaction 
4.1 Pion yields 

The two-component picture of elementary hadronic collisions developed in the preceding 
section can now be used to construct the superposition of target and projectile components also 
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for p+C interactions. The third component given by intranuclear cascading will be shown in 
Sect. |6lbelow to not contribute in the studied xp range. 

Several straight-forward assumptions have to be made in order to quantify this superpo- 
sition: 

- Both the target and the projectile components are assumed to follow the same feed-over 
mechanism as extracted for the elementary collisions as far as relative shape and xp range 
are concerned, as shown in Fig.|9l 

- In terms of pion yields the target contribution has however a weight (z/) where u is the 
number of projectile collisions with the participating target nucleons. The underlying as- 
sumption is here that subsequent projectile collisions lead to the same density distribution 
per participating nucleon. 

- In addition and due to the isoscalar nature of the Carbon nucleus, the density distribution 
for both 7r+ and tt^ is built up as: 



dn \ , 1 / (in V , , , 1 / dn ^^"^ 



, , (7r±) = - -— (7r+) + -{^] (vr-) (8) 

which results as necessary in a tt+Ztt^ ratio of unity over the full xp range. 
- The projectile component on the other hand carries the full imprint of (v) collisions and 
is as such not predictable. Its model-independent quantification and extraction is the main 
aim of this section. In a first instance, it is assumed to be identical to the elementary p-i-p 
collision, 

^ ^-^) = ( (-^) (9) 



with the aim to extract its real shape by inspecting the deviations of the data from this 
simple parametrization. 

With these assumptions, a prediction for the expected pion density distribution in p-i-C 
collisions can be obtained as follows: 

= (-.) + {p-T M (10) 

pred \dXFj,,,^,, \dXpJp,,^ 



For ease of comparison with the data, see Sect. |2l the prediction is referred to the p-i-p 
interaction by forming the ratio 

{dn/dxp)p,^d{xp) 
R,Uxp)- ^^/^^PP(^^) (11) 

This ratio is plotted in Fig. [T21 for three different values of (u) and compared to the 
measurements as a function of Xi?. 

Several points are noteworthy in this comparison: 

- A fair description of the data is obtained in the backward hemisphere at xp towards -0.1 
for the value (z/) = 1.6. The sensitivity of the comparison to (u) is clearly evident from 
the Fig.d 

- This holds simultaneously for the average pion yield and for the 7r+ and vr^ yields sepa- 
rately, both in size and in local slope 

- This is a strong indication for the validity of the two-component picture in this region 
where the target contribution alone should determine the total yields for xp = -0.1 and 
below. 
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- Moving from this xp value towards the forward region, the measurements exceed -Rpred 
up to a; i? ~ 0.2 in a characteristic fashion. 

- For xp > 0.2 the data fall below the predicted values -Rpred- 

The deviation from -Rpred can be translated into absolute particle densities as shown in 
Fig. [13] as a function of x^. 




Figure 13: Difference between measured ratio and predicted yield for a) (tt), b) tt"*", and c) vr 
for (z/) = 1.6 



The obtained density difference follows the shape of the projectile contributions in p+p 
collisions up to x^? ~ 0.1. This is shown in Fig. [141 where the difference is superimposed to 
the projectile charge density distribution for (vr), see Fig. [5] This is a further strong indication 
that the two-component superposition picture produces indeed a valid description of the data. 
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Xp 

Figure 14: Shape comparison of the x p distribution of A (dn/ dxp) in p+C (dashed line) and the 
projectile component {dn/ dxp^m] in p+p collisions (full line) 

The region of reduced particle density at > 0.2 contributes as expected only a very small 
negative contribution to the absolute increase of pion density. Integrating the density differences 
of Fig.[T2]it emerges that the projectile pion yields are enhanced by about 10% with respect to 
p+p collisions 0. This observation carries interesting consequences for the interpretation of 
the corresponding increase of pion production observed in nucleus-nucleus interactions ifTSl . 

4.2 Charge ratios 

The investigation of the tt+Ztt" ratio offers an additional and very sensitive tool for the 
scrutiny of the target-projectile superposition mechanism. As already shown in Sect.|2l Fig.lH 
the overall behaviour of the charge ratio complies with the expectation from a two-component 
picture in the sense that unity is approached in the region of prevailing target contribution, 
xp < -0.1, and the tt+Ztt^ ratio observed in p-i-p collisions is approached in the projectile frag- 
mentation region, xp > 0. 1 . There are, however, detailed deviations from this first-order expec- 
tation. The approach to unity at negative a^i? is not smooth but resembles the structure observed 
in (7r)-i-p collisions dXxp > 0, i.e. in the corresponding isoscalar projectile fragmentation region 
of the elementary hadronic reaction. The tx^Itx^ ratio of p+p collisions is not quite reached in 
the forward region below xp = 0.3. These deviations correspond to a modification of the projec- 
tile charge ratio in multiple collision processes both in the forward and backward regions. A first 
indication for the backward hemisphere is contained in Fig. [13] where the density difference for 
7r+ reaches further out towards negative xp than the one for vr^. This resembles the behaviour 
already found for the elementary collisions. Sect. |3l 

Using the quantitative superposition scheme developed in the preceding Sect. 14.11 for 
(z/) = 1.6, which explicitly contains the tt^/tt^ ratio of unity for the target fragmentation part, 
the deviation of the projectile charge ratio from the one measured in inclusive p+p collisions 
may be extracted. It is shown in Fig.[T5las the ratio {t^^ /'n'^)pfoi/ (^^/^ )inci 

Apparently there is again an increase of the charge ratio in the projectile feed-over region. 
This increase is stronger than the one observed in elementary collisions, Fig.lTOl also shown for 
comparison. It is followed by a decrease of about 4% between xp ~ 0.1 and Xi? ~ 0.25 and 
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Figure 15: Charge ratio of the projectile component in p+C divided by the inclusive ratio in p+p 
collisions 

again by an increase at xp beyond 0.3. This represents an interesting feature of the multiple 
collision process in proton-nucleus collisions in the sense that the projectile charge is correlated 
to both extremes of the fragmentation region. Due to the small contribution in absolute yield 
from these areas, the decrease of the charge ratio in the intermediate region is imposed by total 
charge conservation in view of the increase of projectile multiplicity. 

5 Nuclear aspects 

Up to this point the discussion was performed without recourse to any details concerning 
the manifestation of the Carbon nucleus. Only its isospin symmetry and the general possibility 
of multiple projectile collisions were invoked. It is now necessary to come to a more quantitative 
assessment of a number of nuclear parameters which are essential for the further extension of 
the argumentation. The decisive quantity for most of the nuclear parameters is the nuclear den- 
sity distribution, p(r), which determines the rms radius and the probability of single projectile 
collision, P(l). The second quantity of principle interest is the projectile-nucleon total inelastic 
cross section, or rather its evolution with subsequent collisions. Only by making very distinct 
assumptions about this quantity further important parameters like the probability distribution of 
multiple collisions, -P(i^) and its mean value (u) , already introduced above, may be obtained. 
The quantitative results on several parameters discussed below come from a detailed Monte 
Carlo calculation 1.1 6.1 using the density profile as an input. 

5.1 Nuclear density distribution and percentage of single collisions 

The principle experimental information available for the Carbon nucleus comes from 
electron scattering experiments 1 17 1 which measure the nuclear charge distribution. This is not 
identical with the density profile p{r) as the spatial distribution of charge in the proton has to be 
unfolded. An additional problem resides in the possible asymmetry between proton and neutron 
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distributions. Altogether this leads to a substantial uncertainty about the density profile. Fig.fT6l 
gives p{r) for six different parametrizations: 

1. Saxon- Woods [18|. 

2. Single Gaussian 1 161. 

3. Fourier-Bessel L19I . 

4. Sum of Gaussians |20l. 

5. Fourier-Bessel unfolded for proton charge distribution [ 16.1. 

6. Sum of Gaussians unfolded for proton charge distribution lfT6l . 




Figure 16: a) Nuclear density profile parametrizations; b) rms radius; c) nuclear density at r = 0; 
d) probability for single collision 

In consequence the rms radius, the nuclear density at the origin and the probability of sin- 
gle collisions (which depends on the elementary proton-i-nucleon inelastic cross section) show 
substantial spreads according to the parametrization chosen. This is quantified in the Figs.fT^- 
d for the rms radius, the central nucleon density p(r=0), and the probability of single collision 
P(l). Excluding the first two options which are known to disagree with the experimental evi- 
dence, one may conclude on an rms radius between 2.3 and 2.5 fm and on a fraction of single 
collisions between 0.56 and 0.60. The latter quantity recalls a basic weakness of minimum 
bias proton-nucleus experiments given by the fact that a large fraction of events corresponds to 
"trivial" proton-nucleon collisions. This fraction descends only slowly with increasing atomic 
number A f\V\, as shown in Fig.fTTl 

The contribution of single collisions can in principle be taken out of the experimental 
results since the percentage P(l) is rather precisely determined and the corresponding target 
and projectile components can be safely predicted from elementary collisions. This will be 
demonstrated in the subsequent sections of this paper. An experimentally more efficient way is 
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Figure 17: Probability for single collision and mean number of collisions as a function of 
atomic number 

given by triggering on the presence of "grey" protons f2\\ already at the level of data taking 
which has the additional advantage of permitting the control of centrality. Unfortunately, this 
method could not be applied to the present data due to external constraints on the time available 
for data taking i3j|. 

5.2 The probability distribution of multiple hits and its relation to the total inelastic 
cross section 

Any further specification of nuclear parameters in relation to the interaction of a hadronic 
projectile with the nucleus has to rely on assumptions on the elementary hadron-nucleon cross 
section for multiple collisions. The problem may be visualized in a simple geometrical way by 
looking at a random superposition of A circular disks, each with a surface S of a^P = 31.8 mb, 
as shown in Fig.fTHl 

The surface area covered by one layer of disks is denoted by S 1 . S2 is the area covered by 
two layers, S3 by three layers etc. Scanning the plane with a random distribution of projectiles, 
the mean hit number of layers will be 



The total area covered by the disks is S*! + 5*2 + 5*3 + . . . + S"^. It represents the total 
inelastic cross section a"^^. The total surface of the disks, which is Act^p, is given by 



gi + 2^2 + 3^3 + • • • + ASa 

•S*! + 5*2 + >S'3 + . . . + Sa 



(12) 



AaPP = Si + 2^2 + 353 + . . . + ASa. 



(13) 



The mean number of collisions is therefore defined by the formula 



AcjPP 



(14) 
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Figure 18: Schematic view of the superposition of nucleon-nucleon cross sections in a nuclear 
environment 

This formula is clearly only valid if the disk surface describing the inelastic projectile- 
nucleon cross section is the same in each of the successive layers corresponding to the multiple 
projectile collisions. Any experimental verification of the mean number of collisions provides 
therefore an important constraint on the evolution of the cross section with u. 




Figure 19: a) Probability of u collisions for different nuclear density parametrizations; b) in- 
elastic cross section. The dashed line indicates the NA49 result with the error margin (shaded 
area); c) mean number of collisions; d) mean number of collisions for z/ > 2 
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Always under the above assumption, the distribution P(z^) of the probability to produce 
V collisions in one event, has been determined, using the different nuclear density profiles dis- 
cussed above, via a Monte Carlo calculation lfT6l . The result is shown in Fig. [I9l which also 
gives the integrated quantities (y\^^^, (z/) and the mean number of collisions excluding single 
hits, (y)v>2- 

In principle the measured value of the inelastic cross section at 226.3 mb ± 2% (Fig.figb) 
gives a constraint on the nuclear profile and favours the charge-unfolded distributions (5) and 
(6). This would then correspond to a predicted mean number of collisions of about 1.68. The 
value of 1 .6 extracted from the pion data in the preceding section is indeed close but significantly 
below this prediction as can be seen from the strong dependence on {y) of the density ratios 
shown in Fig.fT2l 

6 Intranuclear cascading 

In the preceding discussion only the principle components of target and projectile frag- 
mentation have been treated. The third contribution to the pion yields due to the propagation and 
interaction of the participant nucleons and of their fragmentation products inside the nucleus, 
will now be investigated. The proper Lorentz frame for the description of this intranuclear cas- 
cading process is the rest system of the target as the cascading products will have typically small 
momenta in this system. The kinematic situation is indicated in Fig. |20l Here lines at constant 
total momentum piab and constant polar angles 9iab for pions in the target rest frame are given 




Figure 20: Kinematic situation in the backward region. With lines are indicated the constant 
total momentum piab (0.2, 0.6 and 1 GeV/c) and constant polar angle Giab (70°, 90°, 118°, 137°, 
160°). The shaded areas represent the coverage of [ 23 1 (far backward region) and NA49 (central 
and forward region). 
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in the xp/pr plane. Pion production is centered at low ptot at ~ -0.15, as expected from the 
mass ratio mjuip. Large values of piab are needed to approach, at forward production angles, 
the region close to Xi? = 0. 

The NA49 acceptance is indicated as a hatched area in Fig. |20l It appears that this ac- 
ceptance matches in its backward low px part the line at 0.6 GeV/c total momentum. At this 
momentum the intranuclear pion production yield is already reduced and only small if any 
cascading contributions are to be expected from this kinematic fact alone. This argument can 
however be quantified by making use of the measurements of Nikiforov et al. 11231 at Fermilab 
with a beam momentum of 400 GeV/c. This important and unique experiment covers the back- 
ward region in the second hatched area indicated in Fig.|20|by measuring identified pion yields 
at constant ©lab between 70 and 160 degrees, leaving only a small uncovered zone with respect 
to the NA49 acceptance. A set of 32 cross sections is available at a;^ > -1.0 for the assessment 
of pion yields in p+C interactions. 

As the experimental layout of the experiment ll23l does not allow a discrimination of off- 
vertex tracks it must however be expected that the complete feed-down from weak decays of 
strange particles is contained in the data at these low total momenta. A feed-down correction 
has therefore been calculated and applied following the methodology developed for the NA49 
data [3]. 

In order to put the Nikiforov et al. data ll23l l into perspective and to allow for a quantita- 
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tive comparison, the given cross sections are related to the isospin averaged prediction for the 
target fragmentation developed and justified in the preceding sections. This is possible since no 
contribution from projectile hadronization can be expected in the covered area. The ratio 

'''''' - {E<Pa/dp%c ^ ^ 

is plotted in Fig. |^ for the 5 measured angles as a function of pt- 

-Rcasc increases strongly with p^ for all angles. This is due to the corresponding increase 
of \xp\ along the lines of constant 6iab together with the increasing dominance of cascading 
products in the far backward region. In fact cascading products are found at < -1 .0 which is 
not accessible for target fragmentation. i?casc therefore tends to diverge towards this limit. The 
arrows given in Fig. ETl indicate the pt value at which xp of -0.85 is reached. 

The dependence on pt at fixed 9iab can be transformed into a dependence on at fixed 
Pt by using eye-ball fits through the data points of Fig.|2T]for interpolation. The result is shown 
in Fig. 122 
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Figure 22: i?casc obtained from ll23ll as a function ofxp 



Within the sizeable error bars of the data points and the interpolation used, a universal 
curve appears which clearly shows convergence towards i?casc = 1 at Xi? = 0. In fact the constraint 
of -Rcasc to 1 at a; = is justified by the fact that relatively large values of piab are requested to 
approach xf = 0, see Fig.|20l As the ratio i?casc gives directly the excess of particle yield due to 
cascading with respect to target fragmentation, it may be argued that the cascading contribution 
is only on a few percent level at xp = -0.1, that is, at the NA49 acceptance limit. The fact 
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that this acceptance limit applies to pt = 0.4 GeV/c or pi^b = 0.6 GeV/c further reduces the 
possible contribution to the NA49 data. The extrapolation into the unmeasured region used in the 
preceding sections, which is oriented towards target fragmentation alone, is therefore justified. 
It should be mentioned here that in the region xp < -0.1 and pt < 0.2, effects concerning 
the tt^/tt ratio must be present due to the final-state Coulomb interaction with the spectator 
protons in the nucleus ll24l . These effects are however expected to be small for the light Carbon 
nucleus. 

In Sect, m above a quantitative prediction from target and projectile fragmentation alone 
has been given extending into the far backward region xp < -0.1. This prediction must now 
be completed taking account of the cascading contribution as quantified in Fig. |22l above. This 
correction has been performed for the pion yield ratio -Rpred shown in Fig. [121 for the charge 
averaged pion yields where it developed a plateau in this region. The complete prediction is 
presented in Fig.l23l 




Figure 23: -Rpred as a function of xp, giving the data points together with the target contribution 
at {v) = 1.6 (full line) and adding the intranuclear cascading part (broken line) 

In this plot, all three principle contributions to the total pion yield, i.e. projectile and 
target fragmentation, and nuclear cascading, are clearly discernible and quantified. The resulting 
overall ratio to the elementary collision exhibits a smooth increase from the far forward to the 
far backward region in xp. This is in agreement with the pseudo-rapidity r] ratio 



Riv) 



{dn/dr]y^{r]) 



(16) 



{dn/dr])P^{r]) 

observed in emulsion work ll25l at 200 and 800 GeV/c beam momentum. Also here a quasi- 
linear increase of R{ri) is observed, from values below 1 at forward rapidity to values in excess 
of 2 in the far backward hemisphere. 
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Contribution from multiple collisions 



As discussed in Sect. |5] and due to the important tail of the nuclear density profile of the 
Carbon nucleus, a large fraction of all events correspond to a single proton-nucleon collision 
in a minimum bias p+C experiment. This probability P(l) is located between 56 and 60% 
corresponding to a reasonable choice of nuclear density distributions (see Sect. |5l and Fig.fT6t. 
Although these single collisions can be accompanied by some amount of nuclear cascading due 
to the single recoil nucleon, this component will not influence the measurements in the NA49 
acceptance region (see Sect.|6l). Given P(l), it is therefore possible to extract from the data the 
content of multiple collisions with u > I which will then rather correspond to a mean value 
(z/) -2.5 (Fig.IIl. 

In a first example the extraction is applied to the charge averaged density ratio R{xf) as 
defined in Eql^labove. In this case the multiple collision component is defined as 



(dn/dxpYrlki^F) R{xf) - -P(l) 



((in/(ixir)P'-(xF) 
and is presented in Fig.|24|as a function of x^. 



(17) 
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Figure 24: Prediction for Rmuh{xF) compared with measurement for minimum bias and 



In comparison to the minimum bias measurement, Rmuitixp) steepens up and reaches the 
values 2.5 at a; = -0.1 and 0.65 at xp = 0.5. It is interesting to regard the limit xp ^ I. Even 
if the suppression of pion density in the far forward region of p+A collisions did not as yet 
find a satisfying explanation, most of the possibilities will finally empty the large x p region of 
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pions from multiple collisions. This is true for projectile energy loss where the scale of xp will 
change, as well as for baryon transfer processes ("stopping"), should baryonic resonance decay 
be the source of forward pions. In these cases Rmuitixp) will approach zero and R{xp 1) will 
reach P(l). This can finally provide a direct measurement of this probability. Fig. l24l contains 
in this sense a tentative extrapolation of i?muit(a;F) towards zero to illustrate the situation which 
will be further constrained by the upcoming publication of proton data in p+C interactions and 
of comparable data from p+Pb collisions. 

A sharp test of this extraction procedure is given by the data sample with centrality con- 
straint also available from NA49 [3|. Here the additional measurement of slow ("grey") protons 
in the lab momentum range below 1 .2 GeV/c provides an event sample where single collisions 
should be rather effectively suppressed. The data points shown in Fig. |24|come from events 
with at least two grey protons detected. As the pion densities for this relatively small subsample 
have been determined in [3 |, the ratio 



{dn/dxpy^^yixp) 
{dn/dxpyv{xp) 



(18) 



can be obtained. Evidently this measurement approaches the predicted behaviour for mul- 
tiple collisions rather closely and proves both the effectiveness of this way of experimental 
centrality control 1 21. . 26.1 and the reliability of the prediction of P(l) from nuclear parameters. 

As a second example for the extraction of improved information on multiple collisions 
the evolution of the mean transverse momentum with multiple collisions will be discussed. The 
situation is shown in Fig. l25l where {px) is shown for p+p collisions 121, minimum bias p+C 
interactions ||3J and for the event sample with at least two grey protons discussed above. 
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Figure 25: Prediction for {pt) in multiple collisions compared with measurement in p-i-p, p-i-C 
minimum bias, and p+C with selected number of grey protons Ug^ey > 2 



Evidently the increase of (pt) from elementary to minimum bias nuclear reactions is 
accentuated imposing centrality control. The multiple collision part can also be extracted from 
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the single collision probability -P(l) via the formula 



where r = {dn/dxpy^ / {dn/dxp)^ = 1/R{xf)- 

The result is shown as the upper line in Fig. 1^51 Again both the measurement with grey 
proton constraint and the prediction from nuclear parameters are coming close. 

A third example of this kind will be discussed in connection with the increase of pion 
yield with transverse momentum in Sect.[Hlbelow. 



Part II 

Pt dependence 



Extending the discussion from px integrated information to pt dependent quantities opens up 
a new and more complex subsurface of the multidimensional phase space as it is projected into 
the double inclusive cross section 

d^a 

f{xF,PT) = E{xf,Pt) ■ -r^{xF,PT) (20) 

dp-^ 

Indeed already the superficial inspection of this quantity in its xf and pt dependence 
for p+p |I3 and p+C 1 3 [interactions reveals significant local structure which has a different xf 
and Pt dependence for the two reactions. In the following sections only a restricted number 
of phenomena can be touched upon, with the aim at bringing out as clearly as possible some 
important features which are specific to p+A collisions and therefore also relevant for the subse- 
quent work on p+Pb and Pb+Pb interactions. The discussion will be conducted again in relation 
to the precise elementary reference given by the NA49 data on p+p collisions O. 

8 Cross section ratio with respect to elementary collisions as a function of pt 

8.1 Definition of the cross section ratio 

The presence of the Carbon nucleus and of the corresponding two-component hadroniza- 
tion mechanism discussed in the preceding sections implies that the straight-forward cross sec- 
tion ratio 

is probably not the best choice for describing the full xf scale. There are in fact two 
important and well quantified overall phenomena which should and can be taken out in the 
comparison: 

- The target pile-up corresponding to the 1.6 participant nucleons which introduces an xf- 
dependent overall upward shift of R in the region a; < 0.1 as elaborated in Sect.|4| 

- The isoscalar nature of the Carbon nucleus which averages both pion charges in the target 
contribution, as opposed to the projectile fragmentation. 
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In order to characterize the projectile-connected pr dependencies in an optimal way, a 
definition of the cross section ratio that takes direct reference to the two-component superposi- 
tion, has been chosen: 

Rpj,{XF,PT) = -? 7 T (22) 

J two-comp\X Fi Pt) 

where /two-comp (a;F,Pr) is defined as: 



/two-comp (a;F,PT) = tixp) (^-/^+ (^F, Pr) + ^f^-^^P^PT)j +p{^f)F'^{xf,Pt), (23) 

with P'^{xf-,Pt) denoting the double differential pion cross section in p-i-p interactions. The 
functions t(xi?) and p{xf) describe the relative contribution of the target and projectile frag- 
mentation to the overall pion yield in the two-component picture. Sects. |3l and lU above, as 
recalled in Fig. |26l 




Evidently, I^xf) ensures that the target contribution amounts to (z/) times the average 
pion yield at < -0.1 and vanishes at = 0.1. On the other hand p{xf) determines the total 
yield dX xf > 0.1 and vanishes towards xf = -0.1. The ratio R^^ therefore contains the salient 
features of the superposition of elementary components in the p-i-C interaction and addresses 
specifically the modifications due to the projectile fragmentation. This is shown in Fig.lTTIfor 
three values of xf and for both 7r+ and vr^. 

Here the full lines represent Rp^ using the xf/Pt interpolation developed for the p-i-p 
and p-i-C |3 1 interactions. As expected, complex p^ and xf dependencies emerge. The pt range 
around the mean values (pr), indicated by arrows in the plots, governs the p^ integrated yields. 
It reproduces the features of this quantity in relation to the elementary collisions as elaborated 
above, namely an increase below xf = 0.2 followed by a zero-crossing and a decrease in the far 
forward region. 
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The relative pt dependence of Rp^ in the vicinity of (pr) shows a rather complicated 
evolution. As the fine structure of /PP(xi?, px) and P'^{xf, Pt) has been connected to resonance 
production and decay [2|, this means that the principle source of pion production changes its 
composition in p+A interactions both as a function of xp and of pt, and in a different way for 
the two pion charges. Here the attention should be drawn specifically to the behaviour at low pr 
governed by the cascading decay of low mass resonances. 

A second region of interest is the evolution towards high values of pr- A significant 



increase of -Rp^ with respect to (px) is evident over the full xp scale as shown in Fig. 1281 which 
provides an overview of Rp^^xp^pT) covering the full range of the measurements. 



8.2 Definition of the pt enhancement 

In order to assess the increase of the pion yields as a function of transverse momentum, 
the measured enhancement factor E^^ is defined as 



T-im 



Rn 



(24) 



This definition relates the behaviour at any p^ to the mean transverse momentum (see 
Sect. |7|). By referring to an integrated quantity which is representative for the pp integrated 
observables discussed in PartHl the detailed structures which are present at low pp are averaged 
out in the reference. The xp dependence of is shown in Fig. |29lfor four values of pp 
extending up to the experimental limit at 1 .8 GeV/c, and using the interpolation shown in Fig.l28l 
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Figure 29: The pp enhancement as a function of xp at several fixed pp values. The shaded 
regions mark the error margins connected to the interpolation scheme |l2l|3|| 

The shaded regions in Fig.|29lcorrespond to the estimation of the statistical error margins 
connected with the interpolation procedure. The enhancement factor increases, from values 
close to unity but significantly different for 7r+ and vr^ at low pp, up to about 2 dX pp = 1.8 
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GeV/c. The increase shows a strong xp dependence, rising from the acceptance limit at xp = - 
0.1 through xf = towards the forward direction. At a;^ > 0.2, 0.3, a decrease of is visible. 

These effects will be discussed in more detail for the high px limit in the following sec- 
tion. 

9 High pt phenomena 

In connection with the discussion of pt dependencies, in particular for proton-nucleus 
interactions, the term "high p^' needs some specification. 

In a first instance, it seems to be a generally accepted notion that hadron production in the 
region of ~ 2 GeV/c is already governed by hard parton-parton scattering. This expectation 
was first developed in connection with ISR data 1 27 1 and is still heavily used in the interpretation 
of very recent RHIC results f2S\. The notion "high" therefore tends to be synonymous with the 
applicability of perturbative QCD, although there is no indication of any serious experimental 
proof for this assumption. 

In a second instance, the early discovery of an "anomalous" px enhancement in proton- 
nucleus collisions, the so-called Cronin effect ||29| . in a corresponding region of transverse 
momentum, seems to point again into the same direction. In fact the only theoretical attempts at 
an understanding of the Cronin effect are based on multiple partonic scattering |l30|. And again 
the present experimental information does not permit in the least to affirm this conjecture. 

In the context of this paper the notion "high" p-r describes the upper limit of pt available 
to NA49 in the p-i-C interaction. The authors consider any relation of this notion to perturbative 
effects as an unfounded speculation. 

9.1 Anomalous nuclear enhancement 

The anomalous nuclear enhancement was first quantified via the A-dependence of 
hadronic cross sections, parametrized as 

ffPA _ ^°(pt) (25) 

The dependence of the exponent a on transverse momentum ll29l is recalled in Fig. |30| 

for 7r+ and tt^. 

The exponent a{pT) rises steeply from values between 0.8 and 0.9 at low pr to a maxi- 
mum of about 1.15 at ~ 4 GeV/c. Above this px the exponent decreases again. The scarce 
measurements and the sizeable error bars at 200 GeV/c beam momentum, i.e. in the immediate 
vicinity of the NA49 data set, do not exclude a sizeable s-dependence. It is also visible in Fig.l30l 
that alpha is slightly larger for vr" than for 7r+. As regards the upper px limit of the NA49 data 
at 1.8 GeV/c, Fig. 1201 the corresponding a^px) is about half-way up from the low-p^ value to 
the saturation limit. It happens to be very close to unity, a value that was believed to be charac- 
teristic of hard parton scattering |29|. The detailed study of the nuclear enhancement in this px 
region gives therefore already valuable information concerning the understanding of the Cronin 
effect. 

The cross sections measured by for Be, Ti and W nuclei may be interpolated to 
Carbon and compared directly to the NA49 data. This comparison is shown in Fig.|2I]for (tt). 

Unfortunately there is only one overlapping data point at 200 GeV/c beam momentum 
in the pr range of NA49, and the consistency of the few data points at 400 GeV/c is not really 
expected due to the strong non-scaling behaviour of the s-dependence in their px range. In this 
context it must be recalled that the data of have been obtained at a constant lab angle of 
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77 mrad which leads to an s-dependent spread of the corresponding xp values between -0.05 
and +0.15, see Fig.|30l In view of the strong xp dependencies of both the target pile-up, Fig.l23l 
and the nuclear enhancement. Fig. |3T1 this leads to systematic effects which are non-negligible 
at the level of precision of the present study. Including however the 20% normalization uncer- 
tainties quoted by [ 291 and the sizeable statistical errors of these measurements, the consistency 
of the data may still be considered as satisfactory. 

9.2 Two basic contributions to the high px enhancement 

Turning back to the NA49 data, it may be claimed that they permit the first systematic 
study of the x p dependence of the Cronin effect, as the data of ||29| only cover a small, s and px 
dependent x p range. The new data at forward rapidity available from d-i-Au collisions at RHIC 
energy suffer from missing acceptance at low px- They will be commented on in a subsequent 
section of this paper. 

Inspection of Fig. l29l shows that the central region of particle production occupies by no 
means a special place, as the enhancement rises smoothly through xp = and has its maximum 
in the region around xp = 0.3 which is far forward with respect to the longitudinal evolution of 
pion yields. This behaviour indicates two important facts: 

- E'^ converges towards unity in the backward region, xp < -0.1. 

- E^^ tends as well towards unity in the far forward region, xp > 0.4. 

In a two-component superposition scheme, the first phenomenon is indeed expected in the 
negative xp limit where the projectile component vanishes with respect to the target contribu- 
tion. In this picture there is no reason to believe that the fragmentation of the participant target 
nucleons shows anomalous behaviour at high px- In other words the target enhancement factor 
^tm-get should be equal to one over the full xp range. The superposition should be characterized 




Figure 30: The exponent a as a function of pp- The dashed line corresponds to pp limit of the 
NA49 data 
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Figure 31: Comparison of the cross section for (vr) as a function of pt at xp = measured by 
||29J and interpolated to Carbon with p+C results from NA49 



by target and projectile components with a specific xp dependence, as presented in Fig.l32b. 

The second phenomenon is of different origin. It is connected to the large probability of 
single collisions, P(l), in minimum bias proton-nucleus interactions, as discussed in Sects. |3] 
and|7]above. It has been shown there that the contribution from single collisions grows as Xi? = 1 
is approached since the pion yield from multiple collisions is progressively suppressed at large 
Xp. As is unity for single collisions, the enhancement should tend towards one for large xp. 
The influence of P(l) on the variation of the measured enhancement with xp can be described 
by a multiplicative weight function e(xp) as shown in Fig. l32b. This function approaches the 
value of 1.6 ~ for 1 and reduces to ~ 0.9 in the backward region due to the yield 

increase from the projectile fragmentation. 



9.3 Extraction of the projectile component 

Concerning the target-projectile superposition, it can be stated that the range of the cor- 
responding target-projectile feed-over is larger than at low transverse momentum (Sect. ^ and 
reaches x^? ~ ±0.3. This indicates an important pt dependence. This is quantified by writing 
down the composition of the measured enhancement from the target (E^^^^^ = 1) and the 
projectile (E^^°^) components: 



rpm 



tixp)E';^s^' + rproj 



/T'proj 

PT 



(26) 



t{xp) + rproj 

where Tproj = {dn/dxp)p^^^{{pT))/{dn/dxpy^{{pT))- Hore t{xp) and p{xf) specify the shape 
of the relative target and projectile contributions as a function of xp. This shape is constrained 
by the boundary conditions p{-0.3) = t(0.3) = and p{xp > 0.3) = l,t(xp < -0.3) = (u), and 
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Figure 32: a) The enhancement factor as a function of x p, b) target and projectile contribu- 
tions as a function of xp Pt = 1-8 GeV/c, and c) efficiency function e multiplying the relative 
contribution from single collisions as a function of a;^ 

the symmetry conditions p(0) = 0.5 and t(0) = {i^)/2 as indicated in Fig. |32| Evidently these 
functions are more similar to the proton than to the low-pr pion feed-over shown in Fig. [TT] 
The enhancement factor for the projectile component can be extracted from Eq.l^as 

= {E^'A^M + rproj) - tixp)) (27) 
The xp dependence of E^^°^ is presented in Fig.l33l 

As an important consequence of this decomposition it appears that the measured nuclear 
enhancement at Xi? = is not fully representative of the phenomenon. It is very substantially in- 
creased if proper account is taken of the "inert" target contribution in p-i-A reactions. This is very 
important for the understanding of the nuclear enhancement in the symmetric A-i-A interactions 
where the target and projectile contribution are equally contributing. From this consideration 
alone it must be expected that the measured nuclear enhancement in A-i-A collisions must be 
larger than in p-i-A collisions. If however properly extracting the projectile component the en- 
hancement should tend to be equal in both types of reaction at equal (z/) per projectile (p-i-A) or 
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Figure 33: The enhancement factors E^, -Ep"j, Ep^°'^>2 ^ function of xj? at p-p = 1.8 GeV/c 



participant (A+A), unless new physics phenomena should be present in A+A interactions. 

9.4 Enhancement due to multiple collisions 

The influence of single collisions on the measured nuclear enhancement is particularly 
important in the p+C case discussed here, as the probability P(l) is about 0.6 and the cor- 
responding enhancement should be unity for this contribution. Taking account of the weight 
function e{xF) introduced above, the enhancement due to multiple collisions can be extracted 
as 



= '^rL.:rr'' (28) 



EP ;°j(x^)-P(l)e(x j.) 
l-Pil)e{xF) 

The application of Eq. |2H] results in a further, very considerable increase of the nuclear 
enhancement over the full xp range as evident from Fig.|33l The final result for E'^^"^^^^ is very 
different indeed from the measured quantity E^. In this context the application of straight- 
forward parametrizations of the A-dependence, as they are applied in the description of the 
nuclear enhancement. Sect. 19.11 above, is to be regarded with great caution. This because the 
probability P(l) depends strongly on A, see Fig.[T7] By relating the effect in p-i-Be to the one 
in p-i-W for instance, one compares cases where P(l) differs by a more than a factor of 3. 
This means that the exponent a^pr) will decrease appreciably once the trivial effect of single 
collisions is taken out. 

As far as A-i-A collisions are concerned, the correct way of comparison would be to 
relate minimum bias p+A to minimum bias A-i-A collisions where the distributions P(z^) are 
equal lfT6l . For centrality selected collisions, this condition is more difficult to fulfill as the 
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probability P(l) stays always greater than zero for A+A interactions even for impact parameter 
zero, whereas it vanishes for central p+A interactions lfT6l . 



9.5 A look at rapidity dependence 



The interpretation of the nuclear enhancement as a px dependent phenomenon necessi- 
tates the use of kinematic variables that are orthogonal in the transverse and longitudinal compo- 
nents, as is the case for the xp/pr combination used here. By describing the same phenomenon 
in the variables rapidity y and pr, a completely different picture will emerge as the use of y 
combines both the transverse and the longitudinal evolution of the particle densities, except, of 
course, at y = 0. This is exemplified in Fig. l34l where the y dependence of the enhancement 
factor 



Ey{y,PT) 



(29) 



Ry{y, {pt)) 

with Ry{y,pT) = {dn/dxFy^{y,PT)/idn/dxFy^{y^PT), is presented at the same four values 
of transverse momentum as shown in Fig. 
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Figure 34: The pr enhancement as a function of y at several fixed px values 

The strong increase of the nuclear enhancement with Xi? is replaced by a flat y- 
dependence up to ?/ ~ 1, followed by a strong decrease for y > 1.5. From this plot alone 
one might erroneously conclude that the Cronin effect vanishes in forward direction, whereas 
is reality one is looking at the longitudinal decrease of particle density in p+A collisions, see 
Sect, m above. 

Experimentally it is often difficult, if not impossible, to fulfill the orthogonality condition 
mentioned above. This is true for all experiments being limited by a polar angle (or rapidity) 
cutoff in forward direction, where the low pr region may not be accessible at all above a certain 
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limit of xf- This is typical of collider experiments. For the BRAHMS detector at RHIC, for in- 
stance, the rapidity limit of 3.2 units corresponds to a px cutoff at 0.8 GeV/c already at a;ir = 0. 1 
where the measured nuclear enhancement is still increasing, and at p-r = 2 GeV/c for xf = 0.25 
where it reaches its maximum. 

10 The trace of resonance decay 

The production and decay of resonances constitutes a much neglected but important field 
of studies concerning soft hadronization. As already pointed out in [2| many details of the 
inclusive pion distributions directly show the imprint of resonance decay. In the context of the 
present paper several of the extracted features concerning the two-component nature of the 
hadronization process will be discussed here in relation to resonance decay. This concerns in 
particular the following experimental findings: 

- The shape and feed-over of the extracted, pr integrated projectile component (Sects. |3l 
and©. 

- Its dependence on particle mass ('Sect. l33t . 

- Its important pt dependence up to the limit of the NA49 data at 1 .8 GeV/c (Sect. |9l). 

As a complete treatment of the manifestations of resonance phenomena exceeds the scope 
of this paper and has to be left to a coming publication, only some basic mechanisms which are 
of straight-forward kinematic origin will be treated here, using a single resonance as a tool for 
explanation. The resonance chosen is the A++(1232) whose inclusive xp and px distributions 
have been determined in an adequate fashion for the present purpose ll3Tl . Its measured, pt 
integrated xp distribution is shown in Fig.|35ltogether with the target and projectile components 
constructed from it, using the baryonic feed-over shape as discussed in Sect. 13.51 The two-body 
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Figure 35: Two-component picture for 
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decay A++ — > p + 7r+ performed for the projectile component of Fig. l35l yields px integrated 
xp distributions of the decay products as shown in Fig.l^ 

A first observation concerns the overall shape of the 7r+ distribution in comparison to 
the extracted inclusive projectile contributions shown in Figs. [5] and [JH Although there are of 
course deviations in detail, as to be expected from the single, low-mass resonance studied here, 
the overall agreement in shape and x^-range is very evident. This concerns specifically the 
limited range of the pion feed-over into the target hemisphere. 

A second observation concerns the difference of this feed-over between proton and pion. 
The mass dependence quantified in Fig. lib is reproduced insofar as the baryon distribution 
reaches much further into the opposite hemisphere than for the pion. 

This effect is of course a straight-forward consequence of the Lorentz transformation 
between the resonance rest frame and the overall cms. The energy factor of this transformation 
reads 



E = ^ml, + g2 (30) 
where msec is the decay particle mass and q is given by 



At the nominal mass of the A++ this corresponds to 0.229 GeV/c which is much smaller 
than the proton mass. The proton therefore receives, in contrast to the pion, a major fraction of 
the resonance momentum in the cms of the interaction. 
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A third observation concerning px dependence is obtained by studying the A++ decay 
for pions in bins of transverse momentum, as shown in Fig.lTTl 

For the same bins used in Fig.|29labove, the distributions are normalized to unity at Xi? = 
in order to calibrate the strong pt dependence of the cross section. Evidently there is a steady 
increase of the range of feed-over with increasing transverse momentum. The shape dependence 
of the Pt integrated xp distributions on the lower integration limit demonstrated in Figs. Eland 
|3] finds thereby a natural explanation. At the upper limit of the NA49 data at 1.8 GeV/c, the 
shape of the pion distribution approaches the one measured for the proton. This effect has been 
used for the extraction of the projectile component of the pr enhancement, Fig.l32l 

This third observation follows again from the resonance decay kinematics. In fact the 
most efficient way to obtain high transverse momentum for low mass secondaries is to ex- 
ploit the Breit-Wigner tail of the resonance mass distribution which reaches, in the case of the 
A(1232), values well above 2 GeV. Correspondingly the decay momentum q (Eq.l?TT) increases 
into the 1 GeV/c region which in turn reduces the mass dependence (Eq. l30t and increases the 
range of the feed-over. It is therefore necessary to take proper account of the resonance mass 
distribution in any quantitative work on hadronization dealing with resonance decay. 

Following these observations, resonance decay as a source of "high p^' production in the 
sense defined above has to be considered very seriously. Using a sum of 1 3 measured resonances 
it has indeed been demonstrated for p-i-p collisions ll32l that the inclusive pion distributions re- 
sulting from two body decays saturate the yields up to = 3 GeV/c and beyond simultaneously 
over the full range of xp. A modification of the resonance spectrum towards higher masses in 
multiple collisions would therefore be sufficient to also describe the anomalous nuclear en- 
hancement in this Pt range, at SPS energies, without invoking partonic or perturbative effects. 
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11 Conclusions 



This paper attempts a new assessment of the p+A interaction with light nuclei at SPS 
energy which is based exclusively on experimental data. The study leads to a number of con- 
clusions which cover a wide range of observables. They may be subsummed as follows: 

• The reference to a set of experimental results obtained in elementary and nuclear inter- 
actions, in particular to the new precision data from the NA49 detector in p+C and p+p 
collisions, allows for a model independent argumentation. 

• This argumentation clearly establishes the composition of the p+A reaction from three 
basic components: projectile fragmentation, target fragmentation, and intranuclear cas- 
cading. 

• These contributions have been isolated and quantified one by one. 

• Intranuclear cascading governs the far backward hemisphere. Its contribution rapidly de- 
creases towards xp = and corresponds to only a couple of percent relative to the partic- 
ipant fragmentation at xp = -OA. 

• For the remaining target and projectile contributions, a clear two component picture has 
been established. The two components overlap in a limited zone around xp = 0. The 
extent of this zone is shown to depend on particle species and transverse momentum. It is 
limited to xp± 0. 1 for the pt integrated pion yield and reaches out to xp± 0.25 at pt of 
1.8GeV/c. 

• Target fragmentation reflects the superposition of participant hadronization, revealing the 

detailed charge dependence introduced by the isospin composition of the nucleus. The 
total pion yield from this source corresponds to an average of 1 .6 participant target nucle- 
ons which is close to but significantly below the number expected from measured nuclear 
parameters. 

• The extraction of detailed features of the projectile fragmentation has been the main aim 
of the study. In relation to elementary collisions a complex pattern of deviations emerges 
which concerns practically all measured quantities. 

• The total pion yield increases by about 10% as compared to the elementary interaction. 
This increase is concentrated near xp = with a shape that matches the density distribu- 
tion of the projectile component. 

• Above Xp ~ 0.2 the differential yield decreases steadily with respect to elementary colli- 
sions as Xp increases, with a specific charge dependence. 

• Using the well defined probability for single projectile collisions, the contribution from 
multiple collisions has been extracted and compared to the centrality controlled p+C data 
subsample. 

• The detailed px dependence of the projectile component has been established in compar- 
ison to elementary collisions. In the neighbourhood of the mean transverse momentum 
the comparison is characterized by important, charge dependent local structure. With in- 
creasing transverse momentum a general increase of pion yield is established. 

• This increase is quantified by an enhancement factor which shows strong xp dependence 
and has been corrected for the inert contributions from target fragmentation and from 
single collisions. 

• The relation to the anomalous nuclear enhancement (Cronin effect) is established and it 
is demonstrated that a new assessment of this phenomenon is mandatory from a purely 
experimental point of view. 

• Finally the importance of resonance production and decay for virtually all studied quan- 
tities is evoked. 
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This study will be continued and deepened by using the data on centrality selected p+Pb 
and Pb+Pb collisions that are also available from the NA49 experiment. 

Acknowledgements 

This work was supported by the Polish State Committee for Scientific Research (1 P03B 
006 30, SPB/CERN/P-03/Dz 446/2002-2004, 2 P03B 04123), the Hungarian Scientific Re- 
search Foundation (T032648, T032293, T043514), the Hungarian National Science Foundation, 
OTKA, (F034707), the PoUsh-German Foundation, the Bulgarian National Science Fund (Ph- 
09/05), the EU FP6 HRM Marie Curie Intra-European Fellowship Program, and the Particle 
Physics and Astronomy Research Council (PPARC) of the United Kingdom. 



41 



References 

[1] H. G. Fischer et al., Nucl. Phys. A715 (2003) 118 
[2] C. Alt et al., Eur. Phys. J. C45 (2006) 343 

[3] Inclusive production of charged pions in p+C collisions at 158 GeV/c beam momentum, 
hep-ex/0606028 (2006) (to be published in Eur. Phys. J.) 

[4] A. Bialas et al, Nucl. Phys. Bill (1976) 461 

[5] D. Barton et al., Phys. Rev. D27 (1983) 2580 

[6] H. G. Fischer et al.. Heavy Ion Phys. 17 (2003) 369 

[7] G. J. Bobbink et al., Nucl. Phys. B204 (1982) 173 

[8] C. M. Bromberg et al., Phys. Rev. D9 (1974) 1864 

[9] T. Kafka et al., Phys. Rev. Lett. 34 (1975) 687 
[10] V. V. Aivazyan et al., Z. Phys. C42 (1989) 533 
[11] S. Uhlig et al., Nucl. Phys. B132 (1978) 15 
[12] K. Guettler et al., Nucl. Phys. B116 (1976) 77 
[13] R. D. Field and R. P Feynman, Nucl. Phys. B136 (1978) 1 

[14] D. Varga, Study of Inclusive and Correlated Particle Production in Elementary Hadronic 

Interactions, PhD. Thesis (2003) Eotvos Lorand University, Budapest 
[15] D. Varga et al.. Heavy Ion Phys. 17 (2003) 387 

[16] A. Rybicki, H.Niewodniczahski Institute of Nuclear Physics, Internal Report No 

1976/PH, http://www.ifj .edu.pl/publ/reports/2006/ 
[17] E. A. J. M. Offermann et al., Phys. Rev. C44 (1991) 1096 

I. Sick, Phys. Lett. 116B (1982) 212 
[18] H. Enge, Introduction to Nuclear Physics, Addison- Wesley World Student Series (1972) 
[19] G. Fricke et al.. Atomic Data and Nuclear Data Tables 60 (1995) 177 
[20] I. Sick, Nucl. Phys. A218 (1974) 509 
[21] M. K. Hegab and J. Hufner, Nucl. Phys. A384 (1982) 353 
[22] A. Bialas et al, Phys. Lett. B51 (1974) 179 
[23] N. A. Nikiforov et al, Phys. Rev. C22 (1980) 700 
[24] J. Burfening et al., Phys. Rev. 75 (1949) 382 

H. Yagoda, Phys. Rev. 85 (1952) 891 

E. M. Friedlander, Phys. Lett. 2 (1962) 38 

N. I. Kostanashvili et al., Sov. J. Nucl. Phys. 13 (1971) 715 
[25] A. Abduzhamilov et al., Phys. Rev. D39 (1989) 86 
[26] S. Afanasiev et al., Nucl. Instrum. Meth. A430 (1999) 210 
[27] M. Delia Negra et al., Nucl. Phys. B127 (1977) 1 

R. P Feynman et al., Phys. Rev. D18 (1978) 3320 
[28] J. J. Gaardh0je et al., Eur. Phys. J. C43 (2005) 287 

B. Z. Kopeliovich et al, Phys. Rev. C72 054606 
[29] D. Antreasyan et al., Phys. Rev. D19 (1979) 764 
[30] J. Kiihn, Phys. Rev. D13 (1976) 2948 

A. Krzywicki et al, Phys. Lett. B85 (1979) 407 

M. Lev and B. Petersson, Z. Phys. C21 (1983) 155 

A. Accardi, hep-ph/0212148 (2003) 

B. Kopeliovich et al., Phys. Rev. Lett. 88 (2002) 232303 
[31] A. Breakstone et al., Z. Phys. C21 (1984) 321 

M. Aguilar-Benitez et al., Z. Phys. C50 (1991) 405 
[32] H. G. Fischer et al., CERN/SPSC 2005-035 



42 



